Jet-quenching in a 3D hydrodynamic medium 
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We study the radiative energy loss of hard partons in a soft medium in the multiple soft scatter- 
ing approximation. The soft medium is described by a 3D hydrodynamical model and we treat the 
averaging over all possible parton paths through the medium without approximation. While the 
nuclear suppression factor Raa does not reflect the high quality of the medium description (except 
in a reduced systematic uncertainty in extracting the quenching power of the medium), the hydro- 
dynamical model also allows to study different centralities and in particular the angular variation of 
Raa with respect to the reaction plane, allowing for a controlled variation of the in-medium path- 
length. We study the angular dependence of Raa for different centralities, discuss the influence of 
hydrodynamical expansion and flow and comment on the comparison with preliminary data. 



Experiments at the Relativistic Heavy Ion Collider 
(RHIC) have established a significant suppression of 
high-px hadrons produced in central A-|-A collisions com- 
pared to those produced in peripheral A-l-A or binary 

\ scaled P+P reactions, indicating a strong nuclear medium 
effect [l|, y| . The origin of this phenomenon, commonly 

[ referred to as jet-quenching, can be understood in the 
following way: during the early pre-equilibrium stage 
of the relativistic heavy-ion collision, scattering of par- 
tons which leads to the formation of dcconfincd quark- 
gluon matter often occurs with large momentum transfers 

\ which leads to the formation of two back-to-back hard 
partons. These traverse the dense medium, losing energy 

\ and finally fragment into hadrons which are observed by 
the experiments. Within the framework of perturbative 
QCD, the leading process of energy loss of a fast par- 
ton is gluon radiation induced by elastic collisions of the 
leading parton or the radiatedgluon with color charges 
in the quasi-thermal medium IHl ■ 

Over the past two years, a large amount of jet- 
quenching related experimental data has become avail- 
able including but not limited to the nuclear modification 
factor Raa, the elliptic flow V2 at high pT (as a mea- 
sure of the azimuthal anisotropy of the jet cross section) 
and a whole array of high px hadron-hadron correlations 
Computations of such jet modifications have acquired a 
certain level sophistication regarding the incorporation of 
the partonic processes involved. However, most of these 
calculations have been utilizing over-simplified models for 
the underlying soft medium, e.g. assuming a simple den- 
sity distribution and its variation with time. Even in 
more elaborate setups, most jet quenching calculations 
assume merely a one-dimensional Bjorken expansion. 



The availability of a three-dimensional hydrodynamic 
evolution code Q allows for a much more detailed study 
of jet interactions in a longitudinally and transversely ex- 
panding medium. The variation of the gluon density in 
such a medium is very different from that in a simple 
Bjorken expansion. A previous calculation in this direc- 
tion [1, 0] estimated the effects of 3-D expansion on the 
Raa- However, this approach treated the energy loss of 
jets in a rather simplified and heuristic manner. Here, 
we shall perform a detailed investigation of the modifica- 
tion of jets in a three dimensionally expanding medium 
within a formalism of p^ . 

Relativistic Fluid Dynamics (RFD, see e.g. [l^, 0, 
[Tsj ') is ideally suited for the high-density phase of heavy- 
ion reactions at RHIC, but breaks down in the later, di- 
lute, stages of the reaction when the mean free paths of 
the hadrons become large and flavor degrees of freedom 
are important. The biggest advantage of RFD is that it 
directly incorporates an equation of state as input and 
thus is so far the only dynamical model in which a phase 
transition can explicitly be incorporated. Starting point 
for a RFD calculation is the relativistic hydrodynamic 
equation 

d^.T^" = 0, (1) 

where T^" is the energy momentum tensor which is given 
by 

T^" = {e+p)U^U'' -pg^\ (2) 

Here e, p, U and g^'^ are energy density, pressure, four 
velocity and metric tensor, respectively. The relativis- 
tic hydrodynamic equation Eq. ([1]) is solved numerically 
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using baryon number ns conservation 
d^inB{T,^i)U'') = 0. 



(3) 



as a constraint and closing the resulting set of partial 
differential equations by specifying an equation of state 
(EoS): e = e{p). In the ideal fluid approximation (i.e. ne- 
glecting off-equilibrium effects) and once the initial con- 
ditions for the calculation have been fixed, the EoS is the 
only input to the equations of motion and relates directly 
to properties of the matter under consideration. Ideally, 
either the initial conditions or the EoS should be deter- 
mined beforehand by an ab-initio calculation (e.g. for 
the EoS via a lattice-gauge calculation), in which case 
a fit to the data would allow for the determination of 
the remaining quantity. Our particular RED implemen- 
tation utilizes a Lagrangian mesh and light-cone coordi- 
nates (r, X, y, if) (r = y/t"^ — z'^), in order to optimize the 
model for ultra-relativistic regime of heavy collisions at 
RHIC. 

We assume that hydrodynamic expansion starts at 
To = 0.6 fm. Initial energy density and baryon number 
density are parameterized by 



^ix,y,v) = emaxW^(a;,y;&)iJ(7?), 
nB{x,y,ri) = UB-ma.xW {x,y;h)H{ri), 



(4) 



where b and Cmax ('^Smax) are the impact parame- 
ter and the maximum value of energy density (baryon 
number density), respectively. W{x,y]b) is given by 
a combination of wounded nuclear model and binary 
collision model [l6j and H{ri) is given by Hijii) = 
exp [-(|?7| - ?7o)^/2cr^ • 9{\-q\ - r^o)] • RED has been very 
successful in describing single soft matter properties at 
RHIC , especially collective flow effects and particle spec- 
tra 0, [13 H, . AH parameters of our hydrodynamic 
evolution [31 have been fixed by a fit to the soft sec- 
tor (elliptic fiow, pseudo-rapidity distributions and low- 
Pt single particle spectra) , therefore providing us with a 
fully determined medium evolution for the hard probes 
to propagate through. 

Let us now discuss the treatment of partons propa- 
gating through the medium: our calculation follows the 
BDMPS formalism for radiative energy loss [lO] using 
quenching weights as introduced by Salgado and Wiede- 
mann fld[2li]. The probability density P(xo, yo) for find- 
ing a hard vertex at the transverse position ro = (a;o, yo) 
and impact parameter b is given by the product of the 
nuclear profile functions as 



P{xo,ya) = 



T^(ro+b/2)TA(ro-b/2) 

TAA{h) 



(5) 



where the thickness function is given in terms of 
Woods-Saxon the nuclear density pa{j^,z) as TA{r) ~ 
J dzp(^A{Tc,z). 

If we call the angle between outgoing parton and the 
reaction plane 0, the path of a given parton through the 
medium ^(r) is specified by (ro, ((>) and we can compute 



the energy loss probability P{AE)path for this path. We 
do this by evaluating the line integrals 

/•oo />oo 

uJc{ro,^)= d^mO and {qL){ro,^)= 

Jo Jo 

(6) 

along the path where we assume the relation 

m =K-2- e3/4(^) (7) 

between the local transport coefficient q{(,) (specifying 
the quenching power of the medium) and energy density 
e. Here, lJc is the characteristic gluon frequency, setting 
the scale of the energy loss probability distribution and 
(qL) is a measure of the path-length, weighted by the 
local quenching power. We view the parameter K as a, 
tool to account for the uncertainty in the selection of 
as and possible non-perturbative effects increasing the 
quenching power of the medium (see discussion in [22l|). 

Using the numerical results of we obtain 

P{AE)path for and R = 2uj'^/{qL) as a function of jet 
production vertex and the angle (here i? is a dimension- 
less quantity needed as input for the energy loss proba- 
bility distributions as defined in [i3|)- The energy loss 
probability P{AE)path is derived in the limit of infinite 
parton energy [lol | . In order to account for the finite en- 
ergy of the partons we truncate P{AE) at AE = Ejet and 
addS{A-Ejet)J^ deP{e). This proced ure is known as 
non-reweighting [llf- We point out that the alternative 
concept of reweighting to our understanding systemati- 
cally overestimates P{AE) for AE ^ Ejet and should be 
disregarded. In fact for a dense medium increasing q and 
employing reweighting leads to an increased escape prob- 
ability whereas increasing q and non-reweighting leads to 
the expected decrease in escape probability, see also [13] . 

Erom the energy loss distribution given a single path, 
we can define the averaged energy loss probability distri- 
bution for a given angle (f> as 



dxo / dyoP{xa,ya)P{AE)path 

OO <J —OO 



(8) 



(this is conceptually similar to the angular averaged dis- 
tribution P{AE))t^a introduced in [23] for central colli- 
sions). 

We calculate the momentum spectrum of hard partons 
in leading order perturbative QCD (LO pQCD) (explicit 
expressions are given in [2^ and references therein). The 
medium-modified perturbative production of hadrons at 
angle can then be computed from the expression 



f 



E ^^^^^ '8>{PiAE))^<S^D}ZH{z, 



(9) 



with Dj°;^f^{z, fj,^) the fragmentation function with mo- 
mentum fraction z at scale 2^, and from this we 
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FIG. 1: Raa for central collisions as calculated in three dif- 
ferent models for the medium evolution 0, [25l . [2^ with the 
overall quenching power scale K adjusted to data. 



compute the nuclear modification function Raa vs. 
action plane as 



re- 



RAAijpT.yA) 



dN'XA/dPTdydcj) 
TAA{h)d<TPP /dPrdydcj) ' 



(10) 



In [2^, l23[ it has been shown that Raa for central 
collisions only constrains a scale, but not the detailed 
functional form of (P(Ai?))T^^. In the approach out- 
lined above, this is manifest in the parameter K which 
we adjust to the data in central collisions. We illustrate 
in Fig. [1] that three different dynamical models, a 2D 
hydrodynamical evolution |25| , the 3D hydrodynamical 
evolution outlined above [3] and a parametrized fireball 
evolution give almost equal descriptions of Raa once 
the scale parameter is adjusted, albeit they require differ- 
ent values of K (the chief reason for this being the differ- 
ent longitudinal dynamics). All three dynamical models 
provide a successful description of the bulk properties 
of the medium at RHIC in central collisions [3, [H, [2^ . 
The parameterized evolution is adjusted in such a way 
that it also describes the Hanbury-Brown Twiss (HBT) 
correlation measurements correctly. Given this success- 
ful description of measured observables within the three 
evolution models the ±50% spread in the values of K 
for the different models of the medium can therefore be 
interpreted as a measure for the systematic error inher- 
ent in the tomographic analysis of jet energy-loss via the 
nuclear modification function Raa- 

However, one may gain predictive power in going to 
collisions at finite impact parameter b. The (j) depen- 
dence of Raa for non-central collisions constitutes a sys- 
tematic variation of path-length within a system with 
fixed overall scale. 

Hydrodyamical models as 0] are able to provide the 
best framework for studying collisions at finite impact 
parameter whereas the application of a parameterized 
evolution model as [1^ to non-central collisions encouters 
the difficulty of how to implement eliptic flow appropri- 
ately. In the following we therefore exclusively resort to 



— P.J, = 6 GeV 

— p.j = 10 GeV 



,< 0.35 




O.25L 



0.5 1 
[rad] 



1.5 



FIG. 2: Angular dependence of Raa for b = 7.5 fm for two 
different values of pr- 
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FIG. 3: pt dependence of Raa in plane (solid) and out of 
plane (dashed) emission at different values of impact param- 
eter b. 



the 3D hydrodynamical evolution [7| which provides a 
very successful framework for the description bulk prop- 
erties also for non-central collisions. 

The average path-length is expected to be smaller for a 
parton emitted in plane as compared to one emitted out 
of plane, and hence Raa is expected to be larger at = 
than at = 7r/2 with the difference in Raa between these 
angles increasing with the initial asymmetry (and hence 
b). Using a simple model for the time-evolution of the 
medium and collective flow effects, it has been shown in 
[2?! that the dependence of Raa is quite sensitive to the 
initial gluon density distribution and temporal evolution 
of the medium. 

Utilizing the previously discussed 3-D RFD model 0], 
we study the angular dependence of Raa for two fixed 
values of at b = 7.5 fm in Fig. [2] and show the pr 
dependence of Raa for emission in plane and out of plane 
at three different impact parameters b in Fig. [31 

As expected, Raa grows for more peripheral collisions 
as there is less soft matter produced to induce energy loss. 
Moreover, there is a smooth angular variation of Raa 
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FIG. 4: Ratio of Raa for out of plane vs. in plane emission 
as a function of pr for four different values of the impact 
parameter b. 

observed, reflecting the underlying medium asymmetry. 
The difference between in-plane and out of plane emission 
grows with impact parameter, at b = 2.4 fm there is 
hardly angular variation whereas at 7.5 fm differences 
are of order 20% (see Fig. |4|). 

It is clear from the observation of elliptic flow that the 
pressure of the hydrodynamical fluid tends to remove 
the initial spatial asymmetry from the almond-shaped 
overlap region in non-central collisions. Thus, on general 
grounds, we may expect that in a dynamical model for 
the medium, the difference between in plane and out of 
plane emission is less pronounced than in an estimate us- 
ing a static medium distributed according to Eq. [Sj On 
the other hand, the main fraction of observed hadrons 
arises from vertices close to the medium surface and the 
expected energy loss per unit time reached a peak 
value early in the evolution [23|, thus it is reasonable 
to expect that a large number of partons escapes the 
medium before the spatial asymmetry is completely re- 
moved. We investigate this competition of timescales in 
Fig. O where we make the comparison with a scenario in 
which we keep the medium static at its initial value. 

Not altering the value K = 3.6 for the quenching scale, 
keeping the energy density at its value in the initial state 
vastly overpredicts the quenching power of the medium, 
leading to Raa of the order of 0.25 for b = 7.5 fm. Thus, 
we have to readjust K to account for the (unphysical) 
fact that we keep the medium static. K = 0.65 leads to a 
good description for central collisions, and employing this 
value also at b = 7.5 fm allows a fair assessment of the 
effect of the hydrodynamical evolution. The azimuthal 
dependence of Raa in the BDMPS-formalism without 
considering expansion effects using a fixed energy density 
profile has also been studied in [28|. 

One finds that indeed the difference between in plane 
and out-of plane emission is reduced by some 7% when 
the hydrodynamical evolution is taken into account prop- 
erly. In addition, a small change of the shape of the dis- 
tribution with pt is induced. Thus, the model shows 



FIG. 5: Raa for in plane and out of plane emission as a func- 
tion of Pt at b = 7.5 fm, assuming a static medium given by 
the hydrodynamical initial state, a static medium with read- 
justed quenching power given by the hydrodynamical initial 
state and the full hydrodynamical evolution. 



sensitivity to the size of the spatial asymmetry in the 
distribution of matter and the timescale at which it is re- 
moved. This complements the information found in low 
Pt V2 which reflects the asymmetry in momentum space. 

However, jet energy loss can also couple to collective 
flow [2^, [sol, thus potentially blurring the relation to 
the spatial asymmetry. In order to assess this effect, we 
implement the effect of flow based on an expression based 
on the definition of q in AdS/CFT. We scale 

9' = g(cosh p — sinhpcos a) (11) 

where a is the angle between flow and hard parton tra- 
jectory and p is the flow rapidity [sij . This expression 
has a straightforward interpretation in terms of the den- 
sity of scattering partners seen by the hard parton per 
unit time/length from the cm. frame. The resulting ef- 
fect on Raa at b = 7.5 fm is shown in Fig. [S) Note that 
the scenario with transverse flow requires an increase of 
K by about 30% to describe Raa in central collisions, in 
agreement with [s^] transverse flow indeed decreases the 
quenching power. 

Transverse flow is in general weak during the timescale 
of jet quenching (i.e. during the early evolution), thus 
there is no pronounced difference between the curves. 
While we observe some change in shape, there is no sig- 
nificant change in the asymmetry between in plane and 
out of plane emission. 

Let us now make a rough comparison with the prelim- 
inary PHENIX data [s^l (since systematic errors due to 
the reaction plane resolution are large, we refrain from 
doing a detailed comparison at this point). While the 
measurement appears to be consistent within errors with 
the calculation for pT > 6 GeV, the data show a trend 
towards a greater angular spread between in plane and 
out of plane emission than found in the calculation. If 
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FIG. 6: Raa for in plane and out of plane emission as a 
function of pr for b = 7.5 fm, comparing the result with and 
without explicit dependence on flow. 

this trend is confirmed, it would indicate that parametri- 
cally the path-length dependence of energy loss is larger 
than and would clearly rule out a linear dependence 
on L (as expected e.g. for elastic collisions as the dom- 
inant mechanism for energy loss). While the depen- 
dence can be understood from coherence length phenom- 
ena (see e.g. [l^l), a dependence on a higher power of 
L has so far no theoretical explanation. The apparent 
disagreement between calculation and data for pT below 
6 GeV is yet another hint that fragmentation is not the 
dominant mechanism for the production of hadrons in 
this momentum regime and other mechanisms, such as 
parton recombination [33 | are of increasing importance 
for lower pT- 

In summary, we have discussed jet energy loss within a 



3D hydrodynamical description of the medium. We have 
made no approximation in computing the average of the 
energy loss probability {P{AE))^ over all possible paths 
through the medium (in particular, we have not tried to 
solve the problem by identifying a typical path). 

Since the nuclear suppression factor Raa is not sen- 
sitive to the detailed form of (P(AiJ))0 as long as the 
quenching scale is adjusted correctly, using a hydrody- 
namical medium does not improve the quality of the de- 
scription of Raa P^f se. However, once we fix the quench- 
ing scale to Raa in central collisions, we gain predictive 
power when going to more peripheral collisions. In par- 
ticular, studying the angular dependence of Raa with re- 
spect to the reaction plane gives systematic control over 
the average in-medium path-length. 

As expected, the angular dependence of Raa reflects 
the spatially asymmetric distribution of soft matter. We 
also observe that the fact that the hydrodynamical evolu- 
tion removes the asymmetry results in a reduction of the 
angular spread as compared to a situation in which the 
asymmetry is kept unchanged. We have also gauged the 
potential impact of the flow fleld on the results. The cal- 
culation describes preliminary data within errors, how- 
ever there is a trend that the data show larger angular 
splitting than the calculation, indicating that the energy 
loss may scale parametrically with a larger power of the 
path-length than L^. It remains to be seen if this trend 
can be conflrmed. 
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